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>  An  analytical  evaluation  of  turbulence-induced  flexural  noise 
in  planar  arrays  of  extended  sensors 
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Large-area,  hull-mounted  conformal  sonar  arrays  typically  employ  extended  sensors  that  are 
configured  to  detect  acoustic  signals  by  means  of  thickness  strains  that  are  induced  by  the 
incident  pressure  field.  In  most  cases,  extended  sensors  also  have  an  appreciable  sensitivity  to 
strains  in  the  lateral  dimensions.  Thus  flexure  of  such  a  sensor  would  induce  a  signal  that  would 
not  be  differentiated  from  that  of  a  target.  Hull-mounted  conformal  arrays  are  evolving  toward 
using  lightweight,  flexible  sensors  and  support  structures;  therefore,  flexure-induced  noise  is  an 
ever  present  concern  This  paper  presents  an  analytical  approach  and  a  general  mathematical 
model  for  the  noise  arising  from  flexure  of  the  array  support  plate  coupled  into  the  array  via  the 
lateral  sensitivity  of  the  sensor.  The  excitation  that  drives  the  flexure  is  assumed  to  be  the 
turbulent  boundary  layer  created  by  motion  of  the  platform  through  the  external  fluid  medium. 

An  analytical  expression  is  derived  for  the  equivalent  plane-wave  spectral  density  for  this  noise 
source.  The  result  is  expressed  in  terms  of  the  frequency  response  function  of  the  plate,  the 
wave-number-frequency  spectral  density  of  the  excitation,  and  the  spatial  filtering  characteristics 
of  the  array.  An  application  is  discussed  to  show  that  predictions  can  be  obtained  in  closed  form. 

PACS  numbers;  43.20.Tb,  43.30.Lz,  43.40.Qi,  43.88.Hz 
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INTRODUCTION 

Large-area  conformal  arrays  that  can  be  mounted  to 
the  hull  of  a  ship  offer  unique  tactical  advantages  over 
towed  arrays.  The  performance  of  such  arrays  is  usually 
limited  by  self-noise  and  platform  noise.  In  the  latter  cat¬ 
egory,  the  noise  induced  by  the  turbulent  boundary  layer, 
located  near  the  hull,  is  a  major  concern. 

Boundary  layer  turbulence  produces  a  random  pres¬ 
sure  field  that  will  be  detected  by  the  array  as  a  noise 
source.  This  is  the  so-called  direct  path  for  flow  noise.  This 
path  for  flow  noise  exists  more  or  less  independently  of 
how  the  sensors  are  supported  and  whether  or  not  they  are 
point  sensors  versus  extended  sensors.  Flow  noise  degrades 
the  signal-to-noise  ratio  but  can  be  reduced  by  using  outer 
decoupler  blankets  that  serve  to  attenuate  the  turbulent 
boundary  layer  (TBL)  pressure  field.  The  use  of  extended 
sensors  to  provide  spatial  filtering  of  the  flow  noise  is  also 
an  attractive  way  to  diminish  flow  noise. 

Secondary  sources  for  flow-induced  noise  can  also  be 
significant.  If  the  structural  support  plate  (SSP)  is  rela¬ 
tively  lightweight  and  compliant,  then  the  TBL  can  induce 
flexure  of  the  SSP,  which  then  serves  as  a  secondary  source 
of  noise.  This  noise  can  enter  the  array  via  direct  flexure  of 
the  extended  sensors  or  as  acoustic  noise  radiated  by  the 
SSP.  The  radiated  component  has  been  addressed  by  other 
investigators.'"*  The  former  source,  flexure  induced  into 
the  sensors,  is  the  focus  of  this  paper. 

The  problem  will  be  modeled  as  follows;  the  SSP,  sen¬ 
sor  array,  and  outer  decoupler  (OD)  are  considered  to 
constitute  a  curved,  layered  shell  with  water  on  the  OD 


side  and  a  vacuum  on  the  SSP  side.  A  vacuum  backing  was 
chosen  because  it  is  simple  to  model  and,  in  addition,  it 
represents  a  worst-case  scenario;  i.e.,  the  case  in  which  the 
SSP  is  backed  by  a  pressure-release  baffle.  The  Corcos’ 
model  will  be  used  for  the  TBL  pressure  spectrum  al¬ 
though  the  theoretical  development  is  quite  applicable  for 
any  model  of  the  wall  pressure  spectral  density.  This  base¬ 
line  model  is  illustrated  in  Fig.  1.  The  formalism  to  be 
developed  makes  no  presumptions  about  the  boundary 
conditions  on  the  plate.  Later  we  shall  assume  that  the 
edges  are  simply  supported  in  order  to  illustrate  a  specific 
application.  Numerous  analytical  studies  of  shell  and  plate 
motion  indicate  that  simple  supports,  for  a  large  shell  or 
plate,  usually  give  results  that  accurately  reflect  the  actual 
response  of  a  plate  supported  in  more  complicated  ways. 
The  fluid  loading  on  the  plate  will  be  included  by  using  a 
rather  simple  model  developed  by  Junger  and  Feit.'“  The 
validity  of  this  model  will  be  established  by  comparing  the 
in-water  displacements  so  derived  with  the  more  exact  pre¬ 
dictions  of  Sandman’s  model." 


3 


2 


I.  ARRAY  RESPONSE  TO  PLATE  FLEXURE 

Typically,  the  extended  sensor  array  is  situated  on  or 
very  near  the  SSP  as  shown  in  Fig.  2.  The  sensors  will  be 
strained  laterally  if  the  plate  flexes  in  response  to  an  exter¬ 
nal  excitation;  therefore,  noise  will  be  generated  if  the  sen¬ 
sor  possesses  lateral  sensitivity.  The  TBL  excitation  will 
also  be  detected  by  the  array,  even  if  the  plate  is  rigid;  this 
is  the  so-called  direct  path  for  flow  noise  as  discussed  ear- 


H 
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FIG.  1.  Configuration  of  baseline  problem. 

Her.  In  order  to  analyze  only  the  flexural  contribution,  it 
will  be  assumed  that  the  sensors  do  not  respond  to  the 
direct  component. 

The  flexural  response  of  a  piezoelectric  plate  has  been 
modeled  by  Ricketts.  “  The  relevant  constitutive  equations 
are 


Ti  —  -1- — /131  Z)3 , 

(la) 

Ti  =  C^iS^  +  C^iSi-hiiDi, 

(lb) 

Ei=  —huSi  —  hiiSi+p^iDi, 

(Ic) 

73  =  74=7-5=0, 

(Id) 

where  T,  S,  E,  and  D  represent  stress,  strain,  electric  field 
intensity,  and  electric  displacement,  respectively.  The  ma¬ 
trix  components  Cjj,  )^3,  and  hjj  are  the  elastic,  dielectric, 
and  piezoelectric  material  constants.’^  The  superscripts  in¬ 
dicate  that  the  designated  parameters  are  held  constant. 


According  to  thin-plate  theory,  the  strains  S,  and  Si 
are  related  to  w,  the  displacement  of  the  plate,  as  follows: 

Si^w  8rw 

Sx=—d-^  and  Si— (2) 

dx~  ay 

where  d  is  the  distance  from  the  sensor  midplane  to  the 
neutral  plane.  Here,  S,  and  Si  are  assumed  to  be  constant 
through  the  thickness  of  the  sensor,  and  is  also  constant 
through  the  thickness,  as  can  be  shown  by  applying  Gauss’ 
law  to  a  dielectric.  Consequently,  is  constant  through 
the  thickness  and  we  can  write 

Ei=V/a,  (3) 

where  V  is  the  voltage  between  the  electrodes  and  a  is  the 
thickness  of  the  sensor. 

Assume  that  the  hydrophones  are  electrically  con¬ 
nected  in  parallel,  which  is  equivalent  to  steering  the  array 
to  broadside.  Typically,  piezoelectric  sensors  operate  in  an 
open  circuit  mode;  hence,  the  total  charge  Qr  appearing  on 
the  electrodes  is  zero.  Therefore,  since  corresponds  ‘o 
the  charge  density,  we  can  write 

Qr=0=  I  Didxdy,  (4) 

,=  1  J  JC,  J  y, 

where  the  integrals  are  performed  over  the  lateral  dimen¬ 
sions  of  each  sensor.  Here,  N  denotes  the  total  number  of 
sensors  in  the  array.  Using  Eq.  (Ic),  we  obtain 

I  £'  =  (5) 

The  first  term  is  independent  of  x  and  y\  therefore. 


where  4  and  ly  are  the  lateral  dimensions  of  an  individual 
sensor.  Combining  Eq.  (5)  with  Eq.  (6)  we  can  write 

{hi\S^-\-hiiSi)A{x,y)dx  dy,  (7) 

—  CO 

where  A(x,y)  is  an  array  sensitivity  function  defined  in 
such  a  way  as  to  turn  the  summation  into  a  continuous 
integration.  For  an  unshaded  array  of  unshaded  hydro¬ 
phones  we  can  write 


A{x,y)  = 


1, 

0, 


if  {x,y)  is  on  a  sensor, 
if  {x,y)  is  not  on  a  sensor. 


(8) 


NO  NOISE  CANCELLATION 

FIG.  2.  Mechanism  by  which  flexural  noise  can  be  induced  via  the  lateral 
response  of  the  sensors. 


Equation  (7)  gives  the  flexural  noise  response  (in 
volts)  in  terms  of  the  strain  components  5,  and  Si-  If  the 
external  excitation  function  is  deterministic,  then  there  will 
exist  unique,  well-defined  strains  that  can  be  computed 
with  thin  plate  or  shell  theories.  On  the  other  hand,  if  the 
excitation  is  a  random  pressure  field  (such  as  TBL),  then 
the  displacement  and  strains  must  be  thought  of  as  sto¬ 
chastic  variables,  which  can  be  represented  by  a  probability 
distribution  function.  In  such  cases,  the  voltage  induced 
into  the  array  will  also  be  a  distributed  variable.  Therefore, 
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in  order  to  properly  assess  the  noise  due  to  flexure,  one 
must  account  for  the  statistical  nature  of  the  excitation, 
which  in  our  case  is  the  turbulent  boundary  layer. 

Before  proceeding  with  the  development  of  a  stochas¬ 
tic  model,  it  is  convenient  to  express  the  noise  sensed  by 
the  array  in  terms  of  an  equivalent-plane-wave  pressure 
field.  This  will  allow  direct  comparison  of  flexural  noise 
with  cir.bient  sea  noise  and  other  noise  specifications.  In  a 
free-field  environment  the  sensors  operate  in  a  hydrostatic 
mode.  An  incoming  plane  wave  of  amplitude  P  will  pro¬ 
duce  an  electric  field  across  the  electrodes  of  a  nonflex¬ 
ing  sensor  of  magnitude 

£'3=(g33+^31+^32)^.  (9) 

where  is  the  transverse  piezoelectric  constant,  and  gj, 
and  g32  are  the  lateral  piezoelectric  constants  for  a  piezo¬ 
electric  slab  that  is  poled  through  its  thickness.  Since  £3  is 
constant  through  the  thickness,  the  voltage  across  the  elec¬ 
trodes  is  V=E^a.  Therefore, 

y^{ag^)P,  (10) 

where  g3,=g33-|-g3i  -fg32  is  known  as  the  hydrostatic  g  con¬ 
stant.  Equation  (10)  allows  Eq.  (7)  to  be  expressed  in 
terms  of  an  equivalent-plane-wave  pressure  field  impinging 
on  an  array  in  the  free  field. 


It.  POWER  SPECTRAL  DENSITY  FOR  FLEXURAL 
NOISE 

The  spectral  density  <l>^p(<u)  for  direct  flow  noise  is 
typically  computed  by  the  following  relationship: 

<J>3,p(«)  =  JJ  d^k  Aik,ks,(o)mk,o) T{k,a)P{k,(o), 

(11) 

where  /4(k,kj,<u)  is  the  array  function  steered  to  kj, 
B(k,oj)  is  the  hydrophone  function,  7’(k,<u)  is  the  transfer 
function,  and  /’(k,<u)  is  the  wall  pressure  spectrum.  The 
essential  features  of  the  derivation  of  this  relationship  are 
established  in  references. Equation  (11)  expresses  the 
direct  TBL  noise  level  at  the  sensor,  but  it  does  not  account 
for  flexurally  induced  noise  that  couples  to  the  lateral  sen¬ 
sitivity.  Therefore  Eq.  (11)  is  not  appropriate  for  the  prob¬ 
lem  under  consideration.  Starting  with  the  first  principles 
that  govern  random  vibration  theory,  an  analogous  expres¬ 
sion  for  flexural  noise  can  be  derived.  The  result  will  be 
formally  similar  to  Eq.  (11),  but  the  interpretation  of  the 
component  functions  will  be  quite  different. 


III.  DERIVATION  OF  THE  POWER  SPECTRAL 
DENSITY  FOR  FLEXURAL  NOISE 

The  following  derivation  of  the  power  spectral  density 
employs  the  notation  and  terminology  found  in  Probabilis¬ 
tic  Theory  of  Structural  Dynamics  by  Lin.'* 

By  combining  Eqs.  (7)  and  (10)  we  can  express  the 
equivalent  pressure  as 


X 


'  X 


{hy^Si-irhiiSi)A{x,y)dx  dy. 


(12) 


As  indicated  previously,  the  excitation  field  is  random; 
therefore,  the  array  output  voltage,  and  hence  the 
equivalent-plane-wave  pressure,  must  be  considered  as  ran¬ 
dom  variables.  The  transverse  displacement  w  and  the  lat¬ 
eral  strains  £|  and  Si  are  also  random  variables.  The  power 
spiectral  density  for  the  equivalent-plane-wave  pressure  can 
be  found  by  taking  the  Fourier  transform  of  the  corre¬ 
sponding  correlation  function  Rp^itj')  that  is  defined  as 


Rpp{tX)^E[PU)PU')],  (13) 

where  £I  ]  indicates  the  expected  value.  Because  P  depends 
linearly  on  5,  and  Si,  we  obtain  from  Eq.  (12) 


5,5,  (M, 


-‘rlhi^^hiiRs  sStt,Ur' ,t') 


+  *32^5j52(f-'.'‘'-'')]'4(r)/l(r'). 


(14) 


where  £5,5, .  '^5,5,  >  ^re  the  cross  correlations  on 

strains.  For  example. 


Rs^Siir,t,r',t')=E[S,{T,t)Si{T',t' )  ].  (15) 

Let  h(T,T',t,t')  denote  the  impulse  response  function  of 
the  plate,  which  is  defined  as  the  displacement  of  the  plate 
at  T,t  due  to  an  impulsive  load  given  by 


P,{r,t)=5(r-t’)d{t-t'),  (16) 

that  is  applied  at  r',f'.  The  principle  of  causality  requires 
that  h{r,r',t,t')  be  zero  when  t<t'.  A  general  excitation 
can  be  represented  by  a  superposition  of  impulses;  likewise, 
the  displacement  can  be  represented  as  a  superposition  of 
impulse  responses  when  the  system  is  linear.  In  addition,  h 
will  depend  only  on  the  difference  t—t'.  Thus  we  can  write 

dV  h(r,r',r-r)P(r',r),  (17) 

where  P  is  the  applied  pressure  and  R  is  the  region  occu¬ 
pied  by  the  plate. 

The  in-plane  strains  are,  therefore, 

d^w 

=  —df  /i„(r.r',/— r)P{r',T) 

(18) 


and 
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d^w 


(19) 


^s^s^(r<r',cj)=d~jj  'J'^p(k,w)C„(r,k, 


X  G*,.{r',k,io)d^k, 


(25) 


where 


where  h,^=d^h/dx^  and  hyy=d^h/dy^. 

Using  Eqs.  (18)  and  09),  the  correlation  functions 
for  strains  can  be  written  as'* 


/?  c  c 


(r,r,r',/') 


J*  dr  J'  ^«(r,s. 


/-r) 


XA;„(r',s',/'-T')/{^p(s,r;s',r'),  (20) 

where  is  the  correlation  on  pressure.  Similar  expres¬ 
sions  are  found  for  jRsjSj  ^tttl  Rs^s^-  The  cross  spectral 
densities  are  obtained  from  the  Fourier  transforms  of  the 
correlation  functions.  For  example,  the  cross  spectral  den¬ 
sity  for  5|  is  given  by 

^5,s,(r,co,r',£u')  =  (^)  f?s,5,(r,t;r',r') 

Xexp-''"'-"'''>r/rdr'.  (21) 


Assuming  the  excitation  to  be  weakly  stationary,  the  cor¬ 
relation  function  for  the  wall  pressure  will  depend  only  on 
the  temporal  separation  t—t'.  In  this  case,  it  can  be 
shown'*  that  the  correlation  functions  for  the  responses 
also  depend  only  on  the  temporal  separation.  Subse¬ 
quently,  the  cross  spectral  densities  depend  only  on  a  single 
frequency  parameter.  Thus  for  stationary  excitations,  Eqs. 
(20)  and  (21)  may  be  combined  to  yield 


d>5_Sj(r,r',6>)=£(^  J  J  d^s' d>p^(s,s» 

XHj^(.r,s,co)H^,^,{r',s',(o).  (22) 

Similarly, 


'I'p^(k,w)  =  |^j  j'J<t>pp(s-s',a))e  *''t/(s-s') 

—  X 

and 


G^^(r,k,aj) 


//^^(r,s,«)e''‘"*  dh. 


(26) 


Similar  expressions  are  obtained  for  4>5  and  5, .  The 
function  (j(r,k,rt))  is  called  the  sensitivity  function.  (As 
before,  subscripts  on  G  and  H  denote  partial  derivatives.) 
This  function  represents  the  structural  response  at  point  r 
when  the  excitation  is  harmonic,  having  wave  number  k 
and  frequency  co. 

For  linear  structures,  Lin,'*  Strawderman,*  and  others 
have  shown  that  the  sensitivity  function  may  be  written  as 
a  superposition  of  normal  modes.  Using  the  method  of  Lin 
we  obtain 


Gj^{T,k,o})  = 


m=l  OJC 

0,  for  re/?. 


for  re/?, 

(27) 


where  /„(r)  is  the  normal  mode  m  for  the  plate  in  vacuo, 
and 


5„,(k)=J  J^/„.(r)e*-'d'r.  (28) 

The  factor  H„,((o)  is  called  the  frequency  response  func¬ 
tion  of  mode  m  and  is  defined  as  the  modal  displacement  of 
the  panel  when  the  excitation  is  a  unit  harmonic  pressure 
having  wave  numbers  corresponding  to  mode  m.  The  func¬ 
tion  S  is  commonly  referred  to  as  the  modal  shape  func¬ 
tion. 

From  Eqs.  (25)  and  (27)  we  obtain 


<I> 


S.S- 


dh’  <I>pp(s,s',cu) 


m,n 


dx^  dx'^ 


Hl{co)H„{co) 


X/f„(r,s,u))//^,^,(r',s',£u),  (23) 

where  the  frequency  influence  function  //(r,s,«)  is  defined 
as 


Lf(r,s,<u)s  r  h{r,s,t)e  dt,  (24) 

«/  —  CO 

and  is  the  Fourier  transform  of  Rpp  with  respect  to 
time.  The  subscripts  on  H  denote  partial  derivatives  with 
respect  to  r  or  r'. 

For  spatially  homogeneous  excitations,  ^pp  will  de¬ 
pend  only  on  the  spatial  separation.  In  this  case,  we  can 
replace  the  integrands  of  Eqs.  (22)  and  (23)  by  their  spa¬ 
tial  Fourier  transforms  to  obtain 


00 

X  JJ  <^pp{k,0))S„,{.k)S*„{k)d^k.  (29) 

—  cc 


Similarly, 

'X 

^5,5j(«'.r',ry)='5i^  Z 

m,n 


dx^  dy'^ 


Hl(co)H„(co) 


X 

X  JJ<I>^^(k,a))S„.(k)5j(k)d'k.  (30) 

—  X 


Equation  (29)  also  gives  <*>5^5^  by  replacing  d/dx  with 
d/dy. 
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Equations  (29)  and  (30)  give  the  cross  spectral  den¬ 
sities  for  strains  in  terms  of  the  modal  response  of  the 
structure  and  the  wave-number  frequency  spectrum  of  the 
excitation.  The  power  spectral  density  for  flexural  noise 
can  now  be  found  by  combining  the  Fourier  transform  of 
Eq.  (14)  with  Eqs.  (29)  and  (30).  Again  invoking  tem¬ 
poral  stationarity,  the  time  Fourier  transform  of  Eq.  (14) 
yields 

<X)  00 

=  (s^)  Jl 

—  oc  —  oc 

-+•  2/»  3 1 ( r,r' ,<u ) 

+  A32‘I>v,(r.'-'.<u)]^(r)^(r').  (31) 

Combining  Eqs.  (29)-(31)  we  obtain 


noise.  The  term  in  Eq.  (35)  is  the  analog  of  the  term 
P(\i,co)  in  Eq.  (11).  It  accounts  for  the  spectrum  of  the 
excitation  field  through  the  relationship  given  by  Eq.  (33), 
where  is  identical  to  /*(k,w).  The  last  term  in  Eq.  (35) 
is  the  analog  of  the  product  of  the  array  and  hydrophone 
functions  in  Eq.  (11).  This  term  depends  on  the  size  and 
location  of  the  sensors  in  the  array.  The  terms  H„(a))  and 
called  frequency  response  functions  (FRF),  are 
not  found  in  Eq.  (11).  These  terms  account  for  the  re¬ 
sponse  of  the  plate  to  each  modal  component  of  pressure  in 
the  spectrum  of  the  excitation  field.  The  FRF  depends 
upon  the  material  properties  of  the  plate  as  well  as  the 
effects  of  water  loading  and  intermodal  coupling. 

The  double  sum  over  modes  in  Eq.  (35)  can  be 
thought  of  as  representing  a  type  of  intermodal  coupling. 
Thus  one  could  have  two  types  of  intermodal  coupling;  ( I ) 
via  the  water  loading  on  the  plate,  and  (2)  via  the  off- 
diagonal  terms  in  Eq.  (35).  These  off-diagonal  terms  arise 
from  cross  correlations  of  normal  mode  strains. 


cJJ  'l>pp(K(o)S„{k)S„ik)d^k 

O^fjr)  aV„(r') 


00  oo 


x/J  a//  A’(a 


31 


dx^  dx'^ 


—  OC  —  CO 


,2  ^fn(r')\ . . 

+  *32 — ^ - q-JT-  U(T)Air  ).  (32) 


If  we  define 


d>p^^(w)=JJ'P^/k.a)5,„(k)</'A:,  (33) 


and 


,  ,,  ,2 

a„„(r,r  >  =  *32 ^+*31*32 

,  ,  ,2 

+  ^31^32  ^  +  *32  . 


then  Eq.  (32)  can  be  written  as 
d 


(34) 


X IJ dvjj aV  a„„(T,r')A(r)A(r'). 

—  00  —  X 

(35) 

Equation  (35)  is  the  central  result  of  this  paper.  It 
provides  the  formal  connection  between  the  wave  number- 
frequency  spectral  density  of  the  excitation  field  and  the 
power  spectral  density  for  the  equivalent-plane-wave  pres¬ 
sure  sensed  by  the  array.  This  relationship  is  the  analog  of 
Eq.  (11),  the  formula  that  is  used  for  direct  path  flow 


IV.  THE  MODAL  FREQUENCY  RESPONSE  FUNCTION 


The  modal  frequency  response  function  can  be 

found  by  solving  for  the  steady-state  motion  of  the  system 
when  a  unit-amplitude,  harmonic  force  is  applied  in  mode 
m.  For  a  rectangular  isotropic  thin  plate,  simply  supported 
and  fluid  loaded  on  one  side,  Lin'*  derives  the  following 
frequency  response  function: 


4 

L^Ly 


Diki+ki)^- 


11(0 


m-(m^,my), 

(36) 


where  L,^  and  Ly  are  the  plate  dimensions,  /r  is  the  mass 
per  unit  area  of  the  plate,  D  is  the  plate  flexural  rigidity,  k^ 
is  the  acoustic  wave  number,  and  p/  is  the  density  of  the 
fluid.  The  modal  wave  numbers  are  given  by 


k„=m^ir/  L^,  = 


and 


(37) 


k„=myTr/Ly,  Wj,=0,l,2,...  . 

Equation  (36)  provides  a  good  model  for  the  fre¬ 
quency  response  function  if  the  mass  and  stiffness  of  the 
isotropic  support  plate  are  much  greater  than  those  of  the 
sensors  and  the  outer  decoupler  and,  in  addition,  the  fluid 
loading  is  adequately  represented  by  the  last  term  in  this 
equation.  When  the  support  plate  is  orthotropic,  and/or 
comparable  in  mass  and  stiffness  to  the  other  components, 
then  the  entire  structure  must  be  thought  of  as  a  composite 
layered  plate  or  shell  as  shown  in  Fig.  3.  In  such  cases,  the 
frequency  response  function  may  be  obtained  using  the 
Donnel  shell  theory  as  generalized  by  Dong’’  to  include 
layered  shells.  If  the  radius  of  curvature  of  the  shell  is 
much  greater  than  the  other  dimensions  of  the  shell,  then 
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FIG.  3.  Cutaway  view  of  the  array  considered  as  a  multilayered 
cylindrical  plate. 

the  Donnel  theory  reduces  to  the  following  equations  for 
the  motion.^® 

Assuming  that  the  plate  is  simply  supported,  the  dis¬ 
placement  can  be  expressed  as  a  summation  of  in  vacuo 
modes;  i.e., 

OQ  CC 

w=  X  X  „  sm{k„  xHin{k„v).  (38) 

m^=I  m^=l  ^  ^ 

Any  externally  applied  pressure  can  be  similarly  written  as 

»  00 

X  X  sin(/:„,  x)sin(*„.  v).  (39) 

m^=l  m^=l  ^ 

The  generalized  applied  force  in  mode  m  is  defined  as 

9/sin(/c„^x)sin(A:„^)£/x’a[v.  (40) 
Using  orthogonality  of  the  in  vacuo  modes,  one  finds  that 


where 

^mnpq  (46) 

^mnpq  (^/v))\Tn(I  L^L^y) .  (47) 

By  examining  Eq.  (44)  we  see  that  each  mode  (m,n)  is 
coupled  to  all  of  the  other  modes  via  the  term 

Leibowitz^^  has  shown  that  the  cross-coupling  terms 
are  much  smaller  than  the  self-impedance  components 


(41), 

The  acoustic  pressure  created  by  the  vibration  of  the  plate 
can  be  expressed  in  like  fashion. 

X  X 

Pa=  X  X  P„  ,n^^m{k„x)im(k„,^).  (42) 

=  1  my.  =  1 

For  a  shallow  layered  shell  it  has  been  shown'”  that 
the  Donnel  theory  yields  ^he  following  equation  of  motion: 

[  Al(^m  +^m  ) +2(  •^12'f  m 

\  y  X  y  X  y 

=  (4/L,Ly)F„^^-P„^„,,  (43) 

where  D^^,  Z)|2,  and  represent  the  flexural  rigidity 
constants  for  a  layered  plate.  These  are  defined'**  in  terms 
of  Young’s  modulus  and  Poisson’s  ratio.  For  a  single  layer 
of  isotropic  material,  these  equations  reduce  to  the  classical 
equation  of  motion  for  a  thin  rectangular  plate. 

V.  WATER-LOADING  EFFECT 

This  problem  has  been  studied  by  many  authors. 

Most  of  the  available  models  have  been  summarized  by 
Leibowitz.^’  The  major  complication  arising  from  water 
loading  is  that  the  orthogonal  in  vacuo  eigenmodes  become 
intercoupled  by  the  water  loading.  That  is,  the  interaction 
of  the  fluid  and  plate  for  one  mode  is  a  function  of  all  the 
other  modes. 

We  shall  adopt  the  Junger  and  Feit  model  as  described 
by  Leibowitz;^’  this  model  was  developed  for  symmetric 
modes  of  the  plate.  However,  the  analytical  result  is  valid 
for  both  symmetric  and  antisymmetric  modes 
(Leibowitz,^^  Table  I). 

The  Junger  and  Feit  model  for  the  water  loading  on 
mode  (m,n)  can  be  expressed  as 

Pmn—^^^^^mnpi^pqi  (44) 

where 


(45) 

I - 

when  k^L^  and  kaLy>'i,  which  are  equivalent  to  k„  L^ 
andk„  Lj,>  1  (the  thin  plate  criteria)  provided  that 
Here  k^  is  the  acoustic  wave  number  cv/c.  This  criteria  is 
satisfied  for  the  SSP  over  the  frequency  band  of  interest. 
Moreover,  Sandman"  has  shown  that  when  moderate 
structural  damping  is  included,  the  cross-coupling  terms 
are  negligible.  For  reassurance  cross-coupling  terms  were 
computed  numerically  and  were  found  to  be  negligible  for 
the  problem  of  interest  to  be  discussed  below. 


^  mnpq 


pcok„k„( 

7? 


(ll 


ryLydr^dy, 

{kl-r],)(k\-ri){k],-r^;)(kl-yi)  ^ 


^^mnpq 


mnpqf 
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TABLE  1.  Water  loading:  Comj>anson  of  approjiimate  and  exact 
methods. 


^mngr 

Frequency  =  8  kHz 
Approximate 
(Ns/m) 

E.xact" 

(Ns/m) 

< 

^1111 

1.368  X  10’ 

(1.357X  10’-yT  .666  »  10') 

'  20  20  20  20 

7.687  X  10’ 

(5.208  X  10’ -y 2.671  A  lO') 

< 

21  21  21 

;5.222xlO’ 

(2.160X  10’-l-)4.%3x  10’) 

/  tn  to  jc'  30 

/T.253X  10’ 

(3.281  X  10’  +  yT.162x  lO’) 

VI.  APPROXIMATE  WATER  LOADING 

When  cross  terms  are  negligible  the  Junger  and  Feit 
model'®  for  the  acoustic  pressure  can  be  written  as 


:■  :i  :! 


'  M  a;; 


;  I,  I'  I  i;  !;  li 
!  :  1  i  :  M  .  f  ;  ; 


i  '  I  ■  !  i  '  '  !  !  I  :  I  / 

'  ■  i  .  I  'I  .  1  I  . 


f'  '  ‘  II  '  i  ;  '  n  I  '  •  'I 

rt'  !•  '1  li  .  !'  1;  '  ' 

1;  O  if  1/  i 


l!  ■ 


i 


X  cooRuiNA  i'i-:  OK  n  Ai  i:  iMi 


ipjio  w„  „ 


■/ti  —  *1 


(48) 


where  py  is  the  density  of  water  and  is  the  acoustic  wave 
number. 

As  a  check  on  the  Junger  and  Feit  approximation, 
some  typical  modes  were  computed  and  the  results  com¬ 
pared  to  the  more  exact  results  of  Sandman. ' '  The  specific 
example  considered  consisted  of  a  0.0254-m-thick  cylindri¬ 
cal  steel  shell  overlaid  with  a  0.0508-m-thick  elastomer 
blanket.  The  lateral  dimensions  of  the  shell  were  L^=4  m 
and  Ly=2  m.  The  radius  of  curvature  of  the  shell  was  5.2 
m.  The  frequency  range  of  interest  was  1  to  10  kHz.  The 
values  used  for  Young’s  modulus,  Poisson  ratio,  and  den¬ 
sity  are  as  follows: 


Elastomer  layer 

£=1.0X10*’  Pa  v=0.49  p=1200  kg/m\ 


Steel  support 

£=2.1  X  10"  Pa  v=0.3  p  =  7900  kg/m'. 


FIG  4.  A  comparison  of  the  plate  displacement  along  the  a  axis  obtained 
using  the  Junger  and  Feit'  approximation  with  the  di.splaeefneni  t>b- 
tained  using  Sandman’s  ’‘metluxl”  for  obtaining  water  loading 


should  it  be  important,  could  also  be  accommodated,  but 
the  numerical  computations  would  become  much  more  dif¬ 
ficult. 

Recall  from  the  structure  of  Eq.  (35)  that  the  double 
summation  over  modes  introduces  an  additional  type  of 
modal  coupling  that  arises  from  cross  correlations  between 
modes.  Thus  one  could  have  two  types  of  modal  coupling: 
(1)  via  the  water  loading  on  the  plate,  and  (2)  via  the 
double  summation  in  Eq.  (35).  We  have  argued  above  that 
the  first  kind  of  modal  coupling  is  negligible  for  our  prob¬ 
lem.  We  shall  see  that  the  second  kind  is  also  negligible. 

Having  established  the  credibility  of  the  Junger  and 
Feit  model  <br  water  loading,  we  are  now  in  a  position  to 
compute  the  frequency  response  function.  Setting 
=  1.  and  using  Eq.  (48)  for  the  water  loading,  we  find  from 
Eq.  (43) 


The  area  density  p  can  be  found  by  taking  the  product  of 
p  and  A.  The  resulting  impedance  values  are  shown  in 
Table  I  for  8  kHz.  The  differences  between  the  approxi¬ 
mate  values  and  the  exact  values  are  found  to  lead  to  errors 
of  no  more  than  1  dB  in  the  final  result. 

The  modal  displacement  as  computed  with  the  radia¬ 
tion  impedance  from  Sandman’s"  numerical  technique 
was  compared  to  a  similar  result  using  Junger  and  Feit’s'® 
approximate  formula,  Eq.  (48).  Figures  4  and  5  compare 
the  results  of  these  two  methods  for  two  representative 
sections  through  the  plate,  one  parallel  to  the  x  axis,  the 
other  parallel  to  the  y  axis.  The  two  methods  give  essen¬ 
tially  the  same  results;  i.e.,  the  differences  are  less  than  the 
resolution  of  the  plots.  These  results  support  the  use  of  the 
simpler  Junger  and  Feit  expression  for  the  purposes  of 
modeling  the  TBL  excitation  of  the  SSP.  The  additional 
developments  that  follow  will  assume  the  Junger  and  Feit 
model  is  being  employed;  however,  these  developments 
could  be  generalized  in  a  straightforward  way  to  include 
Sandman’s  model  for  the  water  loading.  Modal  coupling. 


0.4  - 

t/J 

1- 


n.p 


0  0.5  1  1.5 

Y  COORDINATE  OF  PLATE  (M) 


FIG.  5.  A  comparison  ot  the  plate  displacement  along  the  v  axis  obtained 
using  the  Junger  and  Feit'"  approximation  with  the  displacement  ob¬ 
tained  using  Sandman’s  "method"  for  the  water  loading 
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ipj<o^ 


--77! — 77 


k',  —  k'  —k' 

fft  m. 


(49) 


VII.  THE  ARRAY  FILTER  FUNCTION 

The  in  vacuo  eigenfunctions  for  a  plate  of  dimension 
L  ,  and  Z,„  are 


Ca) 


f„(r)  =:i\n{k^x)sin(k^y). 


(50) 


<b) 


where  m  represents  a  pair  of  modal  indices  that 

take  on  integral  values  from  1  to  x.  From  Eq.  (34)  we 
have 


a„„(r,r')==[A3,ki  + 

+  A32^i,<l/-n(r)A(r') 

=  (^31^m^  +  ^32^mj,n  ^31^:;^  +//52^;;^) 

X/„(r)/„(r').  (51) 


Integration  over  a  single  sensor  located  in  region 
Ri&{Xi,y-x\,y\)  yields 

JJjVjJiVa.. 

22,  22, 


FIG.  b.  Contours  in  the  comptex  plane  used  to  evaluate  the  modal  spec¬ 
tral  density  of  the  turbulent  boundary  layer  pressure  field.  Contour  (a)  is 
used  for  the  diagonal  terms  whereas  contour  (b)  is  appropriate  for  the 
off-diagonal  tern  »»• 


‘  (hy\k\^-\-hiik-„){hy]k\^  +  hi2k‘'„) 

fcos  k„x'-cos  k„  Xi\  /cos  /c„^'  -cos  k„^,' 


(52) 


The  total  array  filter  function  is  obtained  by  summing 
over  all  regions  R,  occupied  by  the  sensors. 


VIII.  THE  MODAL  SPECTRAL  DENSITY  OF  THE 
EXCITATION 

Before  evaluating  <l>„  («)  it  will  first  be  necessary  to 

^  f^mn 

evaluate  the  form  factors  5„(k),  which  are  Fourier  trans¬ 
forms  of  the  mode  shapes;  i.e., 

S„(k)=JJ  e~''‘'''sin(k„^x)sin(k„y)£ix^/y,  (53) 


where  Aj,  denotes  integration  over  the  surface  of  the  com¬ 
posite  SSP.  This  integration  can  easily  be  done  by  parts, 
yielding 


5,„(k)=- 


{ 1  -  exp  I  -  2  { -  m,,7r )  ] } 
k^  -k^ 

k„  ( l-expl-/(kj,L^- ntyTT )  ] } 


X- 


''■my  ''y 


(54) 


The  power  spectral  density  of  the  excitation  <l>^^{k,w) 
will  be  taken  as  the  wave-number-frequency  spectrum  re¬ 
sulting  from  the  Fourier  transform  of  the  wall  pressure 
spectral  model  proposed  by  Corcos.’  Following  Ko  and 
Schloemer,”  this  model  can  be  expressed  as 


P(,Q,(o)(afa2kj) 

(55) 


where  k^=a)/U^  and  is  the  convection  velocity. 

Then  P(0,w)  =04.(1  -|-  y)p]-V*^/o),  where  py  is  the 
water  density  and  is  the  friction  velocity.  The  terms 
a  + ,  y.  a\,  and  02  are  empirically  determined  constants. 
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We  shall  first  evaluate  the  diagonal  terms  of 
=  «).  In  this  case  we  car.  .ombine  Eqs.  (33). 
(54),  and  (55)  to  obtain 

=(i>p^((j)4p^(oj)0(k,),  (56) 

where 


<i>p  (w) 
'  'H 


A;,,  1 1- cost  A, w,77-)jc/A, 

(A,-  A'j  )"jA'-i-  ((f.A  )‘j 

(58) 


and 


(w)  =2 


A;„  [  1  — cos(A,  J  m^Tr)\dky 

(A;-A^  )^[(A,-A,-)-+  (a,A..)-l  ’ 

(57) 


0(A, )  s/’(0.(ij)«|aiA'/  7r. 


t54) 


Consider  Eq  (57),  Tins  integral  can  be  written  as 


4>p  (w)  =  Re  2 


/■ 


I 


a{/.L.  m.rr) 


\dZ 


{Z-k„  )‘(Z+A,„  )-[Z-Av(l+/a,)|(Z^A,(l  -/a,)] 


(MJ) 


where  the  variable  of  the  integration  has  been  replaced  by 
the  complex  number  Z.  This  integral  may  be  evaluated 
using  the  contour  shown  in  Fig.  6(a).  The  numerator  has 
been  replaced  by  the  real  part  of  an  exponential  in  order  to 
avoid  a  pole  where  the  contour  intersects  the  imaginary 
axis.  There  are  simple  poles  at  k„^  and  +  A„  ,  and  a 
second-order  pole  at  Z—kJ^  1  -(-/a]).  Using  the  theorem  of 
residues  we  find  the  following  results: 


<^„  (w)  =  Re 

fm. 


lirki 


1  —exp  i[kpL,(  1  -|-(a| )  — \ 

ai*c-l^c(l+'ai)'-^m,l’  j 


^  / _ 1 _ 

2  \^(,k„^+k^)^  +  (afk,)^ 


(61) 


The  expression  for  (<u)  can  be  evaluated  in  a  similar 

fashion,  the  only  difference  being  that  the  second-order 
pole  will  be  located  on  the  imaginary  axis.  The  result  is 


i  /  .  ,2  l-exp(-a2A,.L^)cosm„n- 

h((o)  =  Ink  - _  r^ia  ,  ,.2 - ^ 


aik^.(aTk^  +  k^  )' 


+ 


-aLy 


(62) 


An  analogous  expression  for  the  off-diagonal  terms 
(m^n)  is 


<l>p  (oj)=4p  {co)<bp  ((o)P(k^),  (63) 

■  /Wff  •  /Tt  *  ^ 

where 


X 


(kl-kl,nkl-kl)[(k,~k,f+(aA)^] 


(64) 


I 

and 

(l_e-  >'''\dk^ 

^  J  *  (  a:  -  A;„')TAf-'  A^TAy+Ta^r-  f’ 

(65) 

These  are  similar  to  the  previous  expressions  except 
that  there  are  now  four  simple  poles  on  the  real  axis  [see 
Fig.  6(b)]  and  the  off-axis  poles  are  now  simple  fHvles.  The 
evaluation  is  straightforward  if  one  writes  the  numerator  as 
I  +  e'”"'  ‘  --  e  ''*>'■  The  first 

three  terms  may  be  integrated  using  the  upper  contour  in 
Fig.  6(b),  whereas  the  fourth  term  may  be  integrated  using 
the  lower  contour.  The  results  will  not  be  presented  be¬ 
cause  they  are  at  least  two  orders  of  magnitude  smaller 
than  the  diagonal  terms  for  the  illustrative  application  dis¬ 
cussed  below.  This  is  to  be  expiected.  for  upon  examination 
of  Eq.  (33)  It  is  seen  that  when  4’^^(k,(j)  is  a  slowly  vary¬ 
ing  function  in  comparison  to  S,„  and  S„,  then  the  contri¬ 
butions  from  the  orthogonal  form  factors  and  S*  tend 
to  cancel  upon  integration.  Thus  the  cross  terms  are  neg¬ 
ligible  except  when  /.,  and  Z.,,  are  so  small  that  and  S„ 
vary  as  slowly  as  Even  if  this  were  not  the  case, 

the  off-diagonal  terms  in  the  array  filter  function  are  small 
for  the  problem  of  interest.  A  numerical  study  of  the  off- 
diagonal  terms  indicates  that  they  are  typically  eight  orders 
of  magnitude  smaller  than  the  diagonal  terms.  The  largest 
contributions  from  off  the  diagonal  are  terms  for  which 
three  of  the  indices  m,,  m,.,  n,,  and  «,  are  equal,  and  the 
fourth  index  is  only  slightly  different  from  the  others.  In 
such  cases,  the  off-diagonal  contribution  is  at  least  two 
orders  of  magnitude  smaller  than  the  diagonal  terms. 

Now  consider  the  product  //,,,( 2u)Z/J(«).  The  diago¬ 
nal  terms  have  peaks  at  frequencies  wher  the  real  part  of 
the  FRF  goes  to  zero.  Figure  7  shows  some  examples  for 
the  case  of  light  damping.  At  the  FRF  resonance  frequen¬ 
cies  the  major  contribution  to  the  modal  sum  comes  from 
only  one  mode.  In  the  case  of  heavy  damping  the  Q  of  these 
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FIG.  7  An  illustration  of  the  diagonal  freituent  y  response  function  for  a 
lightly  damped  support  plate 

resonances  is  greatly  reduced;  thus  more  than  one  mode 
may  contribute.  Even  so,  the  number  of  modes  that  need  to 
be  retained  is  not  large  because  only  nearest  neighbors  con¬ 
tribute  at  that  frequency.  For  example,  as  seen  in  Fig  7,  it 
is  evident  that  to  compute  <I>pp(w)  at  about  1.25  kHz  one 
needs  only  to  retain  modes  that  have  indices  clustered 
about  mode  w,=:/n^.=/i,  =  n,.=  12.  Exclusion  of  the  other 
modes  greatly  reduces  the  computational  complexity  of 
Eq.  (33). 

A  representative  plot  of  off-diagonal  terms 
is  presented  in  Fig.  8.  It  is  seen  that  the 
off-diagonal  terms  may  have  peaks  at  an  FRF  resonance 
frequency,  but  the  amplitudes  of  the  off-diagonal  terms  are 
still  at  least  an  order  of  magnitude  smaller  than  the  diag¬ 
onal  terms. 


FIG  ^  1  he  flexural  noise  spectral  densiiy  as  a  function  of  frequency  for 
the  illustrative  application  discussed  m  the  text.  Results  for  two  types  of 
materials  and  loss  factors  are  shown  The  loss  factors  pertain  to  the 
Young’s  and  shear  moduli  of  the  SSP. 

a„^xdy 

is  at  least  six  orders  of  magnitude  less  for  than  for 

m  =  n.  Therefore,  no  off-diagonal  terms  need  to  be  re¬ 
tained. 

The  following  procedur  was  used  to  evaluate  Eq. 
(33).  First,  the  resonances  of  the  FRF  in  the  frequency 
range  of  interest  were  determined.  For  the  cases  studied 
there  are  about  20  such  resonances.  Then.  <l>^p(w)  was 
computed  for  each  of  these  resonant  frequencies  by  sum¬ 
ming  100  of  the  nearest  neighbor  diagonal  terms. 

X.  AN  ILLUSTRATIVE  EXAMPLE 


IX.  NUMERICAL  PROCEDURES 

Because  the  largest  values  of  the  product 
occur  at  the  resonances  of  the  FRF,  an  up¬ 
per  bound  on  Q{(o)  can  be  found  by  evaluating  Eq.  (35)  at 
those  frequencies  alone.  Only  a  few  modes  about  each  res¬ 
onance  need  to  be  retained  in  the  summation.  Additionally, 
the  off-diagonal  contributions  of  each  term  in  Eq.  (35)  are 
at  least  two  orders  of  magnitude  less  than  the  diagonal 
terms.  Therefore,  the  product 


llHm  Hn*  01 


FREQUENCY  (kHz) 


FIG.  8.  An  illustration  of  off-diagonal  frequency  response  function  for  a 
lightly  damped  support  plate. 


As  indicated  previously,  the  spectral  density  'Pp^ltj) 
represents  the  equivalent-plane-wave  pressure  sensed  by 
the  array.  In  this  form  Qio})  can  be  directly  compared  to 
ambient  sea  noise.  Three  specific  ceises  have  been  studied 
numerically.  All  three  cases  are  dual  layer  laminations  con¬ 
sisting  of  an  SSP  and  an  elastomer  overlayer  of  thickness  of 
0.051  m  (the  OD).  The  array  consisted  of  28x14  pi¬ 
ezorubber  sensors.  The  sensor  dimensions  were  0.095 
X 0.095x0.0032  m.  The  SSP  dimensions  were  £,  =  3.000 
m,  Ly=2.i0  m,  thickness  =  0.0254  m.  The  values  of  the 
empirical  constants  a,,  02,  a  ,  ,  and  y,  that  enter  Eq.  (55) 
are  a  matter  of  some  controversy  in  the  literature.  For  our 
example,  we  shall  use  the  values  suggested  by  Ko  and 
Schloemer,’’  Eq.  (9).  Thus  we  take, 

ai=0.09,  02  =  70,,  0^=0.766,  y=0.389.  (66) 

The  friction  velocity  has  also  been  evaluated  using  the 
following  relationships  suggested  by  Ko  and  Schloemer; 

^j^U,  (67) 

where  U  is  the  platform  speed,  and  Cy,  the  friction  coeffi¬ 
cient,  is  related  to  the  Reynolds  number  /?,  as, 

Cy=0.455(logR,)-"^*.  (68) 

Figure  9  shows  the  results  for  a  0.0254  m  SSP  made  of 
either  steel  or  glass  fiber.  The  platform  speed  is  20  kn.  Two 
values  of  loss  tangent  were  used  for  the  glass  fiber  SSP  to 
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show  the  influence  of  damping.  The  smaller  value  is  more 
representative  of  practical  materials. 

XI.  CONCLUSION 

The  model  presented  above  accounts  for  a  heretofore 
neglected  flexural  no.se  mechanism;  i.e.,  direct  coupling  of 
the  flexure  with  the  lateral  displacement  of  a  planar  sensor. 
The  model  was  developed  for  piezoelectric  sensors  such  as 
rubber-lead  titanate  composites  and  PVDF;  however,  the 
development  could  be  adapted  to  other  sensor  types  such 
as  fiber  optic  sensors. 

An  exact  analytical  expression  for  the  spectral  density 
of  flexural  noise  was  derived  for  the  case  where  the  random 
excitation  field  is  characterized  by  the  modified  Corcos 
model.  In  general,  this  expression  involves  four  summa¬ 
tions,  each  to  infinity,  over  the  normal  modes  of  the  plate. 
However,  it  was  shown  that  only  a  few  modes  need  to  be 
included  for  each  frequency  of  interest.  Therefore,  the 
equation  for  the  spectral  density  is  tractable  for  practical 
applications. 

The  results  indicate  that  a  0.0254-m  steel  SSP  of  light 
damping  and  a  0.191 -m  glass  fiber  SSP  of  heavy  damping 
will  both  perform  satisfactorily.  However,  the  result  for  the 
glass  fiber  SSP  with  light  damping  is  comparable  to  direct 
flow  noise  levels.  It  is  also  seen  that  flexural  noise  increases 
inversely  with  frequency  at  about  the  same  rate  as  sea 
noise.  Therefore,  flexural  noise  arising  from  TBL  excita¬ 
tions  could  be  important  at  lower  frequencies. 

The  model,  as  applied  to  the  example,  did  not  account 
for  the  mechanical  properties  of  the  sensor  nor  for  the 
effects  of  the  inner  decoupler.  The  added  stiffness,  mass, 
and  damping  contributed  by  the  sensor  and  the  ID  will 
lower  ■‘’e  predicted  noise  levels.  Therefore,  the  results  rep¬ 
resent  upper  limits  on  flexural  noise  for  the  various  types  of 
SSPs.  Inclusion  of  the  sensor  in  the  analysis  is  a  straight¬ 
forward  process.  However,  to  include  the  ID,  more  analyt¬ 
ical  development  will  be  required. 

The  dependence  of  the  spectral  density  on  the  array 
design  parameters  can  be  inferred  from  Eqs.  (35),  (49), 
ai.  52)  without  actually  performing  numerical  computa¬ 
tions.  To  minimize  flexural  noise  one  should  do  the  follow¬ 
ing, 

(i)  Minimize  d,  the  distance  from  the  midplane  of  the 
SSP  to  the  midplane  of  the  sensor. 

(ii)  Maximize  and  the  ratios  of  the  hy¬ 

drostatic  sensitivity  to  the  lateral  sensitivities. 

(iii)  Maximize  the  number  of  sensors  in  the  array  and 
the  lateral  dimensions  of  the  SSP. 

(iv)  Minimize  the  gap  between  sensors.  We  see  from 
Eq.  (52)  that  if  the  gaps  were  closed  in  at  least  one  dimen¬ 
sion,  then  there  would  be  essentially  no  noise  due  to  flexure 
of  the  sensors.  However,  there  still  would  be  noise  radiated 
into  the  water  by  the  edges  of  the  SSP. 

(v)  Locate  the  neutral  axis  of  the  flexure  at  the  mid¬ 
plane  of  the  sensor.  For  example,  the  piezorubber  or 
PVDF  type  sensors  could  be  mounted  on  both  sides  of  the 
SSP.  Then,  their  flexural  responses  would  cancel.  How¬ 
ever,  such  a  configuration  would  permit  higher  levels  of 
other  noise  sources. 


The  limiting  assumptions  of  the  model  could  be  ad¬ 
dressed  without  modifying  the  basic  approach.  For  exam¬ 
ple,  the  condition  of  a  pressure-release  backing  to  the  sup- 
p)ort  plate  could  be  replaced  by  a  s[>ecific  impedance 
boundary  condition  which  might  represent  an  inner  decou¬ 
pler  placed  between  the  support  plate  and  the  hull  of  the 
ship.  Another  limiting  assumption  is  that  the  structure  be¬ 
haves  as  a  thin  plate.  This  assumption  manifests  itself  in 
the  model  via  the  relationships  between  lateral  strains  and 
the  normal  displacement  Eq.  (2).  This  limitation  may  be 
removed  by  replacing  Eq.  ( ?.)  with  a  more  exact  expression 
obtained  from  thick  plate  theory. 

In  conclusion,  an  analytical  approach  and  a  mathe¬ 
matical  model  has  been  developed  for  the  case  of  TBL 
noise  induced  via  the  coupling  of  the  lateral  sensitivity  of 
an  extended  sensor  to  the  flexural  vibration  of  the  support 
plate.  This  model  provides  a  means  for  computing  the 
noise  spectral  density  in  terms  of  the  equivalent-plane- 
wave  pressure,  and  thus  provides  a  direct  comparison  with 
ambient  sea  state  noise  levels. 
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